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Abstract —Far-field wireless power transfer (WPT) is a promis¬ 
ing technique to resolve the painstaking power-charging problem 
inherent in various wireless terminals. This paper investigates the 
power transfer efficiency of the WPT segment in future com¬ 
munication systems in support of simultaneous power and data 
transfer, by means of analytically computing the time-average 
output direct current (DC) power at user equipments (UEs). In 
order to investigate the effect of channel variety among UEs on 
the average output DC power, different policies for the scheduling 
of the power transfer among the users are implemented and 
compared in two scenarios: homogeneous, whereby users are 
symmetric and experience similar path loss, and heterogeneous, 
whereby users are asymmetric and exhibit different path losses. 
Specifically, if opportunistic scheduling is performed among N 
symmetric/asymmetric UEs, the power scaling laws are attained 
by using extreme value theory, and reveal that the gain in 
power transfer efficiency is In N if UEs are symmetric whereas 
the gain is N if UEs are asymmetric, compared with that of 
conventional round-robin scheduling. Thus, the channel variety 
among UEs inherent to the wireless environment can be exploited 
by opportunistic scheduling to significantly improve the power 
transfer efficiency when designing future wireless communication 
systems in support of simultaneous power and data transfer. 

Index Terms —Extreme value theory, finite spatial coverage, 
multi-user scheduling, power transfer efficiency, simultaneous 
power and data receiver, wireless power transfer (WPT). 


I. Introduction 
A. Context and Motivation 

C OMPARED to the rapid development of wireless data 
transfer technology, the power transfer technique in 
support of wireless communication systems remained stagnant 
in the past decades. Indeed, wired charging remains the main 
way to feed a device battery, which is a painstaking routine 
for mobile users and can even be a challenge, especially when 
users are on the move or cannot afford interruption in their 
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wireless data service. In particular, when the battery is used out 
and cannot get recharged on time, an abrupt service interrup¬ 
tion does not only degrade the quality of service dramatically, 
but can also yield a loss in users’ data which is unaffordable 
in critical applications. As a result, the shortage of battery 
endurance constitutes a major bottleneck which hinders the 
development of ubiquitous wireless communication systems 
with access to services whenever and wherever needed. 

In order to resolve the power-charging problem inherent 
in current wireless terminals, it is necessary to equip them 
with far-field wireless power charging capability. It is not 
hard to imagine that, if, for example, a mobile phone or a 
laptop can be remotely charged whenever it accesses a wireless 
network, mobile users will never worry about power shortage. 
Clearly, the capability of far-field wireless power charging 
will be infinitely attractive to mobile users and will trigger 
a revolution in the way wireless communication systems of 
the future will function and operate. 

Current techniques of wireless power charging, which are 
mainly based on the principle of electromagnetic inductive 
coupling, are essentially near-field techniques. For instance, in 
the Qi standard, the distance between a charging device, e.g. a 
mobile phone, and its power transfer pad should not exceed 4 
cm (1.6 in) Q. Also, the power transfer pad must be connected 
to an electrical outlet in a wired way. Obviously, compared 
to traditional wired charging, the near-field wireless power 
charging can reduce neither the users’ concern on the shortage 
of battery endurance nor the users’ labour in the charging. 

To extend the distance of wireless power transfer (WPT), 
magnetic resonance and/or radio frequency (RF) based far- 
field power transfer are promising. The concept dates back 
to the late 1800s. Indeed, the earliest prototype of magnetic 
resonance based WPT was demonstrated by Heinrich Hertz 
in 1887 Fig. 24]. In practice, Nikola Tesla attempted to 
wirelessly transmit power by magnetic resonance at Colorado 
Springs, CO, USA, in 1899. In this experiment, however, no 
evidence was collected on whether any significant amount 
of power would be available at any distant point. The first 
successful RF based power transfer was performed by Harrell 
V. Noble at the Westinghouse Faboratory, which was demon¬ 
strated to the general public at the Chicago World’s Fair of 
1933-1934 j^. Eater on, RF based far-field WPT technology 
found extensive applications where no human interaction is 
needed, for instance, aeronautics and aerospace. Among them, 
the most famous experiment was performed by William C. 
Brown in 1963, who wirelessly powered a helicopter fly¬ 
ing at 60 feet above the ground level. Another ambitious 
WPT experiment was the Stationary High Altitude Relay 
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Platform (SHARP), which was performed by Communications 
Research Centre of Canada in 1987. In this experiment, RF 
based WPT was used to provide energy from ground station 
to SHARP at about 70,000 feet altitude. Then it was not until 
2008 that simultaneous wireless transport of power and data 
in communication systems was proposed, by Lav R. Varshney 
at MIT 0 - Although magnetic resonance and/or RF based 
WPT techniques have already found many applications in 
practice 0 0 Chap. 4], potential health risks of human 
exposure to electromagnetic fields constitute a major factor 
which needs further R&D efforts to achieve the full potential 
of these techniques and a widespread commercial deployment. 
This safety issue is beyond the scope of this paper and the 
interested reader is referred to the up-to-date survey in Q. 
Since magnetic resonance based WPT and RF based WPT 
have different design methodologies, in the rest of the paper 
the term WPT is limited to the latter, i.e. WPT via radiowaves, 
unless otherwise stated. 

Figure [2 shows a wireless system in support of simultaneous 
power and data transfer, where the block diagram of a user 
equipment (UE) consisting of a specific power receiver and 
a traditional data receiver is sketched. In this system, a 
conventional base station (BS) exchanges data with the UE 
while a newly deployed power beacon radiates power to the 
UE. In practice, power beacons can be deployed separately 
from BSs or be integrated into BSs. In particular, power 
beacons can be integrated into femtocell BSs, since they both 
serve nearby users within the area of radius on the order of ten 
meters. With such a dedicated power beacon, power charging 
cables are not compulsory anymore and UEs can be remotely 
charged without the users’ intervention. Although this system 
is promising for next-generation wireless communications, 
there are few related papers in the open literature with most 
of them dedicated to calculating the achievable power versus 
data rate region, from an information theory perspective. 
Specifically, the original work 0 focused on non-fading 
scenario, which was extended to frequency-selective fading by 
considering a coupled-inductor circuit problem in Q. Lately, 
the authors of 0 studied how to design power/data receivers 
in a separate or an integrated way. Based on the separate 
structure in optimal switching between the power receiver 
and the data receiver was addressed in Go), to achieve various 
tradeoffs between power transfer and data reception. 

In practice, in order to avoid interference between the power 
receiver and the data receiver of a UE, these units should 
work in different frequency bands. Actually, since far-field 
WPT is still in its infancy, so far no particular frequency band 
has been specified yet for this technique. Thus, so far the 
industrial-scientific-medial (ISM) radio bands remain the only 
option in prototype development. For example, in 0, 5.8 
GHz and 2.45 GHz were used for power transfer and data 
transmission, respectively. Accordingly, at the UE side, the 
power receiver and the data receiver are separately designed, 
as shown in Fig. Generally speaking, systems in support of 
simultaneous power and data transfer at different frequency 
bands are known as out-hand systems whereas those working 
at the same frequencies are known as in-hand systems. Clearly, 
out-band systems are free of interference between the power 
and data transfer segments. 



Fig. I. A wireless system in support of simultaneous power and data transfer. 


B. Efficiency of Ear-Field Wireless Power Transfer 

Figure shows the inner structure of a WPT segment 
in an out-band system (corresponding to the link from the 
power beacon to the power receiver shown in Fig. [T]), which 
consists of five subsystems: (1) DC-to-RF conversion, (2) 
retrodirective beamforming at the transmit antenna array, (3) 
free-space transmission, (4) receive antenna array and (5) 
rectifier. Thus, the end-to-end power transfer efficiency, 7 ^e 2 e, 
can be expressed as the product of 77 ^, V/ G [1,5], where 77 ^ 
denotes the power transfer efficiency at the subsystem. 
That is, 

P ^ 

A out 


'^e2e 




( 1 ) 




where Pin and Pout denote the input and the output DC powers 
of the WPT segment, respectively; and where 771 = Ptx/Pin, 
V2 = -Perp/Ax, V3 = Anc/Arp, V4 - Ax/Anc and 
?75 = Aut/A-x, as shown in Fig. Alternatively, 7]e2e can be 
decomposed as the product of the transmit efficiency (771 x 772 ), 
the free-space propagation efficiency ( 773 ) and the receive 
efficiency (774 x 775 ). 

In the state-of-the-art, in order to improve the end-to-end 
power transfer efficiency, researchers’ efforts were focused on 
enhancing the transmit efficiency, the receive efficiency, or 
both. To obtain higher transmit efficiency, the transmit antenna 
array must be designed such that the side lobes of its beam 
pattern are reduced to the lowest acceptable level and that 
its main lobe keeps spillover losses to a minimum. To this 
end, electronically steered phased arrays with retrodirective 
beamforming have taken the place of traditional horn antennas 
and emerged as the most reliable technique to guarantee 
an accurate beam steering in WPT applications. On the 
other hand, attaining higher receive efficiency was attempted 
through the design of high-performance rectifying antennas 
(i.e. rectennas), which convert the incident RF power back to 
DC. A comprehensive survey on the design of transmit and 
receive antenna arrays dedicated to WPT applications can be 
found in GD It is noteworthy that omnidirectional antennas 
must be avoided in WPT applications in order to improve the 
power transfer efficiency and, in particular, decrease the risks 
of human exposure to electromagnetic fields. 

Despite the efforts to improve the transmit efficiency (771 x 
772 ) and/or the receive efficiency (774 x 775 ), the main deter- 
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Fig. 2. Block diagram of a WPT segment in an out-band system. 


minant of the end-to-end power transfer efficiency remains 
7 / 3 , i.e. the free-space propagation efficiency. In fact, recalling 
the well-known Friis transmission equation, for an effective 
radiated power (ERP) Perp at a power beacon (defined as 
the product of transmit power and antenna gain), the incident 
power Pine at a receive antenna array is exponentially decayed 
with the Euclidean distance between the power beacon and the 
power receiver (cf. Eig. [^. Without loss of generality, setting 
the path-loss exponent to the urban environment value of 4, 
the power transfer efficiency can readily be shown to be much 
smaller than 1 %, given that the distance between the beacon 
and the receiver is on the order of ten meters. Therefore, 
improving the power transfer efficiency is a prerequisite for 
the success of far-field WPT. 

C. Contributions and Organization of the Work 

In this paper, we exploit multi-user scheduling, which 
was proven to be effective in improving data transmission 
efficiency in wireless communication systems, to improve the 
power transfer efficiency in propagation space (i.e. 7/3 shown 
in Pig. 1^ by utilizing the channel variety among the UEs. 
Moreover, since the rectifier is the core component of a power 
receiver, its non-linearity effect is accounted for and the output 
DC power. Pout, shown in Pig.[^ is quantified. Por simplicity. 


the receive antenna array is assumed ideal and, thus, 7/4 is 
normalized to unity. 

To study the effect of channel variety among UEs on the 
efficiency of far-field power transfer, we first look into the 
homogeneous scenario where UEs, described as symmetric in 
this case, have a similar distance from the power beacon, 
and we analyze the time-average output DC power when 
either round-robin scheduling or opportunistic scheduling is 
implemented for the power transfer to the UEs. Then, we 
investigate the heterogeneous scenario with asymmetric UEs, 
which are uniformly distributed in the considered annulus, and 
analyze the average output DC power with the said scheduling 
strategies. The results of this study reveal the effects of the 
inner and outer boundaries of the system coverage, multi-user 
scheduling and channel variety among UEs, on the power 
transfer efficiency of the WPT segment in an out-band system. 
In particular, three major contributions of this paper are: 


1) Implementing round-robin scheduling, the channel va¬ 
riety among asymmetric UEs is demonstrated to be 
equivalent to a distance scaling factor, compared to the 
case with symmetric UEs. Also, the said factor is shown 
to be determined by the size of the inner and outer 
boundaries of the power beacon coverage and by the 
path-loss exponent; 

2) Performing opportunistic scheduling while considering 
N symmetric UEs, the power scaling law is derived and 
reveals that the scheduling gain is InA^, compared to 
the power transfer efficiency pertaining to round-robin 
scheduling; 

3) Applying the same opportunistic scheduling as above 
among asymmetric UEs, the power scaling law is an¬ 
alytically attained and discloses that the scheduling 
gain is as large as N (number of UEs), compared to 
the power transfer efficiency pertaining to round-robin 
scheduling. Therefore, the channel variety among UEs 
can be exploited to improve the far-field power transfer 
efficiency significantly. 


To detail the paper’s contributions, the remainder of the 
manuscript is organized as follows. Section ^ describes the 
system and the channel models. In Section [ffl the inner struc¬ 
ture of the UE power receiver is sketched and the amount of 
instantaneous output DC power at the receiver is formulated. 
Section [Iv| is devoted to the scenario with symmetric UEs 
and analyzes the time-average output DC power and power 
transfer efficiency with different scheduling policies. Then, 
Section |V] focuses on the scenario with asymmetric UEs. 
Simulation results and discussions are presented in Section [Vl| 
Concluding remarks are provided in Section VII followed by 
mathematical tools and derivations relegated to appendices. In 
particular, the third contribution mentioned above was reported 
in the accompanying conference version GD- 

Notation: Eor a complex variable x, operators 5R{x}, ^{x}, 
\x\ and arg{x} denote its real part, imaginary part, amplitude 
and phase, respectively. E{x} indicates the statistical expec¬ 
tation of real random variable x. Constant j = ^/^ denotes 


the imaginary unit. The functions T{x) = t 




dt and 


r(a, x) = ^ dt, Va,x > 0, refer to the Gamma 

function and the complementary incomplete Gamma function. 





















4 


IEEE TRANSACTIONS ON SIGNAL PROCESSING, ACCEPTED EOR PUBLICATION 



Fig. 3. The WPT segment in the simultaneous power and data transfer 
system illustrated in Fig. ^ where the power beacon (PB) in the center of 
the circular area is steering power to the user terminal, UE^, i G [1, A^]. 
All terminals are located in the annulus of inner radius rex and outer radius 

r* net • 


respectively. The Landau notations f{x) = 0 {g{x)) and 
f{x) rsj g(^x) are defined as linia^^oo |/(^)/^(^)| < cxo and 
f{x) /g{x) = 1, respectively. Finally, Pj' stands for 
the time-average output DC power with respect to the joint 
states X and y, where x G {rrs, os} refers to the scheduling 
policy (‘rrs’ for round-robin and ‘os’ for opportunistic) and 
y G {horn, het} to the status of UEs (‘horn’ for the homoge¬ 
neous case and ‘het’ for the heterogeneous case). 


II. System and Channel Models 

As depicted in Fig. we consider the WPT segment in an 
out-band system where the power beacon steers power to the 
user equipment, UE^, Vi G through a directional 

antenna array. The power beacon is located in the center of a 
circular coverage with radius r^et, with an exclusion zone of 
radius rgx around it such that all UEs are in the far-field of the 
power beacon. The exclusion zone can be seen as a prohibited 
zone where users are not allowed for human safety reasons due 
to the high radiated power of power beacons (more details on 
the safety issue will be examined in Remark 4). 

At the power beacon, the radiated RE signal at time slot t 
is expressed as 

s{t) = (2) 

where Pgrp = E{5^(t)} denotes the average effective radiated 
power at the transmit antenna array, / refers to the carrier 
frequency, and x{t) is a complex baseband signal of bandwidth 
B Hz and unit power (i.e. \x{t)\‘^ = 1). Also, it is assumed 
that wireless channels experienced by s{t) keep fiat during 
each transmission slot. 

A basic wireless channel model between the power beacon 
and receiver UE^ consists of large-scale path loss and small- 
scale multi-path fading. Mathematically, the instantaneous 
multiplicative channel gain Gi{t) between the power beacon 
and UE^ at time slot t is given by 

Gi{t) = /3d-'^{t)\hi{t)\^ (3) 


where ^ > 0 is a constant scaling factorQ(ij(t) e [rex, rnet] 
is the distance between the power beacon and UE^, a > 2 
denotes the path-loss exponent, and hi{t) stands for the 
complex channel coefficient. 

III. Power Receiver 

With respect to the receiver, Eig. illustrates a typical 
rectifier-based power receiver, where the incident RE power. 
Pine, is first converted by the rectifier into direct current 
(DC). Then, a boost converter is cascaded to provide the 
required step-up from typical rectenna voltage (usually, tens 
to hundreds millivolts) to typical battery voltage (usually, 2 
to 4 volts). As shown in the dashed block of Eig. the 
rectifier is composed of a Schottky diode cascaded with a 
low-pass filter (LPE), thanks to the low power-loss and the 
fast switching/recovery time of Schottky diode, compared with 
ordinary P-N junction diodes. 

In light of and 0, the received signal at the power 
receiver of UE^ can be shown to be expressed as 

yi(t) = ^2PerpGi{i) 5R {a:(t) } + n,(t), (4) 

where Oi{t) G [0,27r) and n^(t) denote the phase shift 
introduced by the propagation channel and the additive white 
Gaussian noise (AWGN) at the receive antenna of UE^, re¬ 
spectively. Without loss of generality, rii (t) is assumed to have 
zero mean and variance 6‘^. Notice that, for the considered far- 
field WPT segment, the average effective radiated power, Pgrp, 
is generally much larger than thermal noise variance, S‘^. As 
a result, 0 can be well approximated by using 

Viit) « ^2PerpGi(i) 3? {x{t) } , (5) 

Then, for an easier mathematically tractable form of yi{t), 0 


can be reformulated as: 


Viit) 

« \/23? j 




= V2?fi\ 

0PerpGi(i) eJ[<^i(t)+arg(a,(t))]gi2^/t 

} 


= V2 A(t) cos [27r/t + 0(f)] , 

(6) 

where the amplitude A{t) and the phase 0(t) are given by 



(7) 


(pit) = arctan , 

^ ^ yi{t) 

(8) 

with 





yi{t) = 

\rpGi{t) COS [6»i{q + arg(a;(i))], 

(9) 


yoX) - \IPerpGi{t) sin[6»i(i) + arg(a;(i))]. 

(10) 


Since the current-voltage characteristic of a Schottky diode 
can be described by the well-known diode law, in view of the 


Tn general, parameter (3 in Eq. ^ stands for the log-normal shadowing 
effect introduced by objects obstructing the propagation path between the 
power beacon and the UE. This effect influences the local-mean powers at 
the UE . tha t is, short-term averages to remove fluctuations due to multi-path 
fading [T^ Section 2.4], which is accounted for by the parameter defined 
immediatMy after Eq. QZ) of the paper. 
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Fig. 4. Block diagram of a typical rectifier-based power receiver. 


sinusoidal input voltage given by ([^, the output current of the 
Schottky diode in Fig. can be readily given by 

(exp (^) - l) ’ (11) 

where Is denotes the reverse saturation current of the diode, p 
is the quality factor of the diode and Vt refers to the thermal 
voltage. The value of p typically varies from 1 to 2 depending 
upon the fabrication process and semiconductor material and, 
usually, is assumed to be approximately equal to 1. On the 
other hand, Vt = kT /e, where k is the Boltzmann constant, T 
is the working temperature in Kelvin, and e is the magnitude 
of the electronic charge. Subsequently, by using the Taylor 
series expansion of the exponential function in ( pT] ), the output 
current can be rewritten as 

OO J 

= ( 12 ) 


After passing through the LPF immediately after the Schot¬ 
tky diode, the harmonic components at kf (k > 1 ) of i{t) 
shown in © are removed and only the DC component 
remains at the output of the rectifier (cf. FigJ^. Therefore, by 
ignoring the high-order (k > 2 ) terms in and performing 
some trigonometric manipulations, the DC output of the 
rectifier can be readily given by 




(13) 


With the input current idc(^), the converted power to be 
stored in the rechargeable battery is in general linearly pro¬ 
portional to the value of idc{t), with a conversion coefficient 
0 < ^ < 1 (^ > 0.85 in practice 0 ). Thus, the instantaneous 
output DC power Pout-i(^) at UE^ during transmission slot 
t, in the unit of watt or equivalently joule/sec, is given by 




^-^s-^erp 


Gi{t) = cl3dr-{t)\hi{tp, 

(14) 


where c = ^IsPerp/{pVT)‘^• Clearly, ( p^ demonstrates that 
the instantaneous output DC power at UE^ during each trans¬ 
mission slot is proportional to the path loss which 

is the fundamental reason why the WPT efficiency degrades 
significantly with the propagation distance. 


Recalling the five subsystems of Fig. 13 it is noteworthy that 
the mathematical model starting from Q to •El captures the 
physical processes where the signals propagate over space and 
are processed by the rectifier in a power receiver illustrated in 
Fig.|^ Hence, the power transfer efficiency described hereafter 
is equivalent to the product of ps and ps whereas pi, p 2 and 


?74 (cf. Fig. IJl are not accounted for 0 

In order to investigate the effect of the distance di on 
the power transfer efficiency in a multi-user context, in the 
sequel we consider two different scenarios in conjunction 
with different multi-user scheduling strategies for the power 
transfer from the beacon. 

Remark 1 (Rectifier and rectifying theory). The typical 
rectifier-based power receiver shown in Fig. is a normal 
wave detector suitable for relatively weak power receiver 
and/or energy harvester. It enables mathematical tractability 
without any side effect on the major findings of the paper. For 
medium and high power-consumption applications, a single¬ 
shunt full-wave rectifier can be used, which consists of a diode 
and a capacitor connected in parallel to a A^/4 distributed 
line, with Xg being the effective wavelength of an input 
radiowave 0 Fig. 3.67]. For the general rectifying theory 
and the state-of-the-art of various rectennas, the interested 
reader is referred to 

IV. Scenario A: Homogeneous Users 

In this scenario, UEs are assumed to be uniformly dis¬ 
tributed along a circle centered by the power beacon or located 
closely to each other, widely known as clustered users (e.g. 
several tablets on a table while the PB is adsorbed onto the 
ceiling above, or users in a small coffee shop). In such a 
case, UEs have almost the same distance from the power 
beacon and, thus, they get the same average power as per 
•El- Accordingly, for ease of presentation, in this section 
the distance dfit) is abbreviated as d, regardless of the user 
index (i) or the time variable (t). Actually, if we set d = 1, 
this idealized scenario reduces to the case where the path-loss 
effect is ignored, as usually assumed impractically in the open 
literature related to WPT. 

Next, in order to study the effect of path-loss on the 
WPT efficiency and, in particular, provide a performance 
benchmark for later development, two different scheduling 
strategies among symmetric UEs are considered and their 
corresponding average output DC powers are analyzed. 

A. Round-Robin Scheduling 

In round-robin scheduling, UEs take periodic turns to access 
the medium of service. Under this policy, all UEs get an equal 
share of the available resource, such as time, and thus have 
the same performance. 

^For more details on the way to improve the subsystems’ efficiencies ryi, 772 
and ? 74 , the interested reader is referred to and references therein. 
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By recalling that di (t) = d for all UEs in the homogeneous 
case, the instantaneous output DC power at UE^, given by 
( p^ , can be rewritten as 

Pout-i{t) = cl3d~“\hi{t)\'^. (15) 

Assuming Rayleigh model for the effect of multi-path fading, 
the term \hi{t)\‘^ in ( p3] ) has an exponential distribution such 
that the probability density function (PDF) and the cumulative 
density function (CDF) can be expressed as follows: 

f\h,it)\^{x) = Xexp (-^] , (16) 

^\hi{t)\'^{x) = 1 - exp ’ (17) 

where = E{|/i^(t)p}, Vi G denotes the average 

small-scale multi-path gain. By virtue of (p3])-(p^, it is easy 
to derive the PDF and the CDF of Pout-i(^), namely, 

= (18) 

When round-robin scheduling is implemented, each UE gets 
the same average output DC power. In view of this power 
value can be easily computed as 

_ POO 

xfp^^^_.(^t){x)dx = cl3ald-'^, ( 20 ) 

Jo 

which implies that the average output DC power increases 
linearly with larger multi-path gain (cr^) but decreases ex¬ 
ponentially with the transmission distance (d), as expected. 
Although the result in ( [20| ) is straightforward, it will serve as 
a benchmark for the subsequent complex cases. 

B. Opportunistic Scheduling 

Although round-robin scheduling policy ensures absolute 
fairness among UEs, in terms of access to the system re¬ 
sources, it is not very efficient. For instance, for UEs which 
experience deep fading, the amount of instantaneous output 
DC power can be extremely small. In order to enhance the 
power transfer efficiency, opportunistic scheduling, which was 
originally proposed to improve spectral efficiency, can be 
exploited |T^ . In this case, assuming that a UE cannot get 
fully charged when scheduled, the UE with the maximum 
instantaneous output DC power is chosen to be charged during 
each transmission slot. To this end, a transmission slot is 
divided into two sub-slots: one training sub-slot and one 
power-transfer sub-slot. During the training sub-slot, each UE 
sends its received power level back to the power beacon by a 
reliable feedback mechanism, and the power beacon chooses 
the one which has the highest power level to be charged in 
the following power-transfer sub-slot. With such a multi-user 
scheduling strategy, the index of the chosen UE is given by 

I = arg max Pout-i(i), (21) 

i=l, ••• , N 

where Pout-i(^) denotes the instantaneous received power of 
UE^, defined in ( p3] ). By recalling the results of order statistics 


theory, the CDF and the PDF of the maximum output DC 
power Pout -2 (^) of the chosen UE^ is readily shown as 

= ( 22 ) 
= N (23) 

where and are given by ([T8]( and 

respectively. Accordingly, the average output DC power 
can be computed by 

_ pOO 

Phom = d:r. (24) 

Although substituting ( p^ and ( p^ into ( [^ and perform¬ 
ing some mathematical manipulations leads to a closed-form 
expression, the resultant formula is very complex and has to be 
evaluated numerically. In order to offer illuminating insights 
into the output DC power, we instead derive the limiting 
distribution of the maximum output DC power Pout- 2 (f), as 
N ^ oo. It IS noteworthy that the limiting distribution cannot 
be obtained by directly applying ^ oo in ( [22| ) since, for 
any Ppo^^_^^(t)(x) < 1, ( [22| ) reduces to 0 sls N ^ oo and, 
thus, the CDF given by ( [22| ) is a degenerate distribution. 

In the following lemma, the asymptotic theory of extreme 
order statistics is exploited to attain a non-degenerate distribu¬ 
tion for Pout -2 (f) and, then, the power scaling law is derived. 

Lemma 1 (The limiting distribution of the maximum instan¬ 
taneous output DC power among symmetric UEs). If the PDF 
and CDF of the output DC power at UE^, Mi G 
are respectively given by and then as the total 

number of UEs, N, becomes large asymptotically, the limiting 
distribution of the maximum instantaneous output DC power 
Is of Gumbel type, namely, 

lim +6ia;) = exp (-e"®) , (25) 

N^OO 

where the positioning parameter ai and the scaling factor bi 
are given by 

ai = d~^ In (26) 

hi = cP(j‘f^d~^. (27) 

Proof: The proof is provided in Appendix ■ 

Recalling the fact that the mean of the Gumbel function is 
Euler’s constant, i.e. 0.5772 • • • GZl p. 298], with the resulting 
Lemma it is straightforward that the average output DC 
power can be computed as 

phom = d-»(lnN + 0.5772 • • • )■ (28) 

In light of ( |28| l, we immediately have the following theorem. 

Theorem 1 (Power scaling law for symmetric UEs). When 
there are N UEs which have the same distance (d) from the 
power beacon, if the UE with the maximum instantaneous 
output DC power is chosen to be charged during each 
transmission slot, then as N becomes large asymptotically, 
the average output DC power at each UE scales as 

PhgO™ ~ cjdal d~" In N. (29) 

Proof: When the value of N becomes large asymptoti¬ 
cally, ignoring the Euler’s constant in the parentheses of ( |^ 
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yields the desired ( |^ . ■ 

By definition, the power transfer efficiency pertaining to op¬ 
portunistic scheduling is given by while that per¬ 

taining to round-robin scheduling is expressed as P^^°™/Perp. 
Accordingly, the ratio of their power transfer efficiencies can 
be shown as P^g°™/P^^o™. Consequently, by comparing ( [29] ) 
with ( [20] ), we conclude the following corollary. 

Corollary 1 (Power transfer efficiency for symmetric UEs). 
The opportunistic scheduling among N symmetric UEs yields 
a scheduling gain of InA^ in terms of the power transfer 
efficiency, compared with the round-robin scheduling. 

Remark 2 (On the multi-user scheduling in the context of 
power transfer). In the context of power transfer, the received- 
power based multi-user scheduling, given by Eq. ( [21] ), is 
essentially equivalent to the channel state information (CSI) 
based multi-user scheduling in the context of data exchange 
in conventional communication systems. However, unlike the 
data receiver where the CSI can be obtained by the built- 
in channel estimation unit, it is not necessary for the power 
receiver to estimate CSI and, thus, the feedback information 
from the power receiver to its target power beacon is only the 
received power level. On the other hand, due to the relatively 
very short duration of the training sub-slot in comparison with 
the subsequent power-transfer sub-slot, the energy obtained at 
each terminal during the training sub-slot is negligible. 

V. Scenario B: Heterogeneous Users 

When UEs are uniformly distributed as depicted in Eig. [^ 
the distance di{t) between UE^ and the power beacon is a 
random variable with respect to the user index i and time 
variable t. In such a case, the average output DC power of 
different UEs vary dramatically with their respective distances 
from the power beacon; a scenario referred to as asymmetric 
UEs. In this section, taking the distance variety into account, 
we investigate the power transfer efficiency of the system by 
means of the time-average output DC power. Eirst, we start 
with the scheduling according to the round-robin policy. 


A. Round-Robin Scheduling 

When UEs are uniformly distributed in the annulus of inner 
radius rgx and outer radius r^et (cf. Eig.|^, the distance di{t) 
of UE^ from the power beacon is a random variable and its 


CDE and PDE are respectively given by 



to 

1 

if rex < a: < rnet, 

Otherwise; 


Pi{t) (x) 

) ^2 _ ^2 ’ 

\ ' net ' ex 

lo, 

(30) 


f 2x 

if “^ex ^ ^ ^ ^net 5 

Otherwise. 


fd.it) (x) = ■ 

} r 2 _ ^2 ’ 

\ ' net ' ex 

io, 

(31) 

In light of ([T^, ([Tt]), ([30]) and ([ 3 T]), and performing 

some 


mathematical manipulations (details provided in Appendix ]B]), 
the CDE of the instantaneous output DC power given by ( |l^ 
can be found to be expressed as 

FPont-iit) (^) 


= 1 - 


^net 


^net 


[rlet exp (-r 2 x) - rl^ exp (-n a;)] 


cj3a\ 


a 


r 1 H—, ria; - r IH- ,r2X 


a 


(32) 


where n = r^J{cPal) and r 2 = r"^J{c/3al). 


When round-robin scheduling is adopted, all UEs get the 
same average output DC power and, with the resultant ( [3^ , 
it can be computed as 


phet 

rrs 


0 

oo 

2 _ 2 / exp (-r 2 x) - Tgx exp (-ri x)] da: 

^net ^ex J 
0 


^2 _ ^2 
' net ' ex 




X J X ^ 
0 


r [ 1 n - r fi r2 X 


dx 


rlet - ^ex V ^2 


f 


a{cl3al)i 


2c^al 


{a - 2)(r2^, - rlf) 


ri J 




rlx)^ j 

(33) 

/ 2 — 0; 2 —o;\ 

\ (^ex ^net j ? 

(34) 


where Eq. (6.455.1)] was exploited to derive ( [3^ . Thus, 
( [34] ) establishes the relationship between the average output 
DC power and the system parameters in a very simple and 
explicit way. Eurthermore, it is general in practice that rnet ^ 
Tex and, thus, ( [34] ) can be approximated as follows: 


phet 

rrs 


a- 2 V^net; "" 




Ti 


2 

ct-2 



T 2 


2 


(35) 

(36) 


By comparing ( [3^ with ( [^ , it is observed that: 1) The 
term T\ of ( [36] ) corresponds to the maximum average output 
DC power of the symmetric-user case with the shortest dis¬ 
tance Tex from the power beacon; 2) The term T 2 of ( [3^ is a 
power scaling factor (with respect to Ti) that reveals the effect 
of path-loss variety among UEs on the average DC power. 


On the other hand, if we set ( [^ equal to ( [20] ), we may 
find a symmetric-user case, whereby the users’ distance from 
the power beacon is d, and users get the same average output 
DC power (accordingly, the same power transfer efficiency) 
as in the heterogeneous case. Specifically, d is determined by 


cPaf^ d ^ 


c(3alr 


—a 

ex 


2 


r’ex 


2 


5 


Q (-2 


net 


(37) 
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which yields 


d = 


a ■ 


-fey 


T 3 


(38) 


Clearly, the term T 3 in ( [^ is a distance scaling factor (with 
respect to the shortest distance rgx) which reflects the effect 
of path-loss variety among UEs on the average DC power. In 
other words, if round-robin scheduling is adopted, the average 
output DC power of asymmetric UEs is equivalent to that of 
the symmetric ones at a distance rex^s from the power beacon. 


B. Opportunistic Scheduling 

When opportunistic scheduling is applied among the asym¬ 
metric UEs, it yields higher scheduling gain than the above ho¬ 
mogeneous case. To show this, by using again the asymptotic 
theory of extreme order statistics (but completely different 
from the proof of Lemma [^, we obtain the following result. 


Lemma 2 (The limiting distribution of the maximum instanta¬ 
neous output DC power among asymmetric UEs). If the CDF 
of the instantaneous output DC power at UE^, Vi G 
of the heterogeneous scenario is given by ( [32| ), then as N 
becomes large asymptotically, the limiting distribution of the 
maximum instantaneous output DC power Pout- 2 (^) is of 
Frechet distribution, namely, 

^ , (39) 


where the scaling factor 62 is given by 


ho = 


‘^cNpal / 

a-2 feneti “ 


(40) 


Proof: The proof is provided in Appendix ■ 

Since the path-loss exponent a > 2 holds in practice, 
the exponent of x in the right-hand side of ( [^ is no 
less than —1, i.e. — - > —1. In this case, the mean of 
the Erechet distribution approaches infinity, i.e. the integral 


rOO 

Jo 

cannot 


1 — exp 


(■*:*) 


dx does not converge. Therefore, we 

to 


use a similar methodology as in Section IV-B 


attain the power scaling law. To proceed, we derive the upper 
and lower bounds on the natural logarithm of the maximum 
instantaneous output DC power, and then obtain 

the power scaling law. 

Specifically, by virtue of Lemma the CDE of Pont-iit) 
can be explicitly shown as 


Pr |pout-.(f) < (^) re-“4 = exp , 

(41) 

for X > 0 and N -fe 00 . By recalling the fact that the 
logarithm function In x is a monotonically increasing function 
with respect to x > 0 , can be reformulated as 


Pr |lnPout-?(t) < In ^ ^^2 (ru) + 

(-x-S), (42) 


= exp I 


Now, taking x = In A/", it is clear that 
limjv^ooexp (lnA^)“«^ = 1 and, thus, ( |42| ) implies 

lnPout-t(0 < In ( f Ux“ I + In A" + Inin A" 

\ a - 2 \rnet/ I 

(43) 

almost surely when N ^ oo. 

On the other hand, (|4^ can be equivalently rewritten as 


Pr < InPout-z(t) > In 


2cl3al 
a — 2 



+ In iV + In (r 


= 1 — exp 



), (44) 


for X > 0 and N ^ oo. Now, taking x = 1/In A/", it is 
clear that limjv^cxD exp ( 1 /In A^)“ « ^ = 0 and, hence, 
( |44| ) implies 

In Pout-t(f) > In f + In A - In In A 

\a-2 \rnetj J 

(45) 

almost surely when N oc. 

Einally, combining ( |4^ and ( [45] ) yields 



^ 1 


(46) 


almost surely when N ^ oc. Eor completeness of exposition, 
the above result is summarized in the following theorem. 


Theorem 2 (Power scaling law for asymmetric UEs). When 
there are N UEs uniformly distributed in the annulus of 
inner radius rex cind outer radius if the UE with the 
maximum instantaneous output DC power is chosen to be 
charged during each transmission slot, then as N becomes 
large asymptotically, the logarithm of the average output DC 
power scales according to 

(47) 

Next, by comparing © with ( [35] ), it is evident that 

ln/^-ln(AJ^). (48) 

Afterwards, by using a similar approach to attain Corollary 
we have the following corollary. 


Corollary 2 (Power transfer efficiency for asymmetric UEs). 
The opportunistic scheduling among N asymmetric UEs yields 
a scheduling gain of N in terms of the power transfer 
efficiency, compared with the round-robin scheduling. 

Comparing Corollary with Corollary it is clear that the 
path-loss variety among asymmetric UEs leads to higher WPT 
efficiency. Specifically, if opportunistic scheduling is adopted, 
the power transfer efficiency increases from \nN in the ho¬ 
mogeneous case (no path-loss variety among UEs) to N in the 
heterogeneous scenario (with path-loss variety among UEs). 
Hence, by implementing opportunistic scheduling instead of 
the traditional round-robin scheduling at a power beacon, far- 
field power transfer efficiency can be significantly improved. 


for X > 0 and N ^ oo. 
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TABLE I 

Parameter Setting Used in the Simulation Experiments 


Symbol 

Definition 

Value 

Unit 

f' ex 

Radius of Exclusion Zone 

2 

m 

^net 

Radius of Circular Coverage 

30 

m 

Is 

Reverse Saturation Current of Schottky Diode 

1 

mA 

N 

Number of User Terminals 

1-50 


P 

Transmit Power of the Power Beacon 

43-53 

dBm 

V'jp 

Thermal Voltage 

28.85 

mV 

a 

Path-Loss Exponent 

4 


5 

Shadowing Effect 

1 


p 

Quality Factor of Schottky Diode 

1 



Coefficient of Energy Conversion 

0.85 


-1 

Average Multi-Path Gain 

1 



Remark 3 (On the fairness among users in the heterogeneous 
scenario). The coverage area of a power beacon is generally 
very small, compared with that of a BS. For instance, in 
the simulation setting detailed in Table the radius of the 
coverage area of the considered power beacon is set to 30m 
and a typical UE is located 10m away from the beacon, 
whereas the radius of the BS cell is usually as large as 
500m. Therefore, the power levels received from the beacon at 
different UEs are on the same order. As a result, the fairness 
among heterogeneous users in the context of power transfer 
is almost guaranteed. 


VI. Simulation Results and Discussions 

In this section, we present and discuss simulation results, in 
comparison with numerical results pertaining to the previously 
developed analysis. The parameter setting used is summarized 
in Table |I| In particular, the radii of the exclusion zone and the 
circular coverage of the power beacon are set to 2m and 30m, 
respectively. Also, a typical user is set to be 10m away from 
the power beacon, and its achieved power transfer efficiency is 
used to compare the effectiveness of different WPT schemes 
with or without scheduling. In the following, we first discuss 
simulation results obtained using round-robin scheduling. 


Remark 4 (Safety levels of human exposure to RF elec¬ 
tromagnetic fields). According to the IEEE Standard C95.1- 
2005, for safety levels with respect to human exposure to 
RE electromagnetic fields, the permissible exposure level from 
2GHz to lOOGHz in a public environment is lOW/mf [19 


p. 27]. Typical frequencies that are used for prototype devel¬ 
opment of far-fieId WPT systems are 2.45GHz and 5.8GHz 
(because they belong to the ISM bands), fall exactly into 
the aforementioned frequency range. In our simulation setting 
(Table^, if the effective radiated power (Perp) cit the power 
beacon takes the maximum value 53dBm, i.e. 200W, the power 
densities at the distances from the power beacon d = 2m (the 
closest allowable), 10m (distance of a typical UE) and 30m 
(edge of beacon coverage), are 3.9SWjw?, 0.159Wjw? and 
0.0177WIw?, respectively, by recalling that power density is 
computed via Pq = Perp/(47r(i^). Evidently, these densities 
are smaller than lOW/rnf and, thus, permissible in practice. 


A. Round-Robin Scheduling 

Figure illustrates the average output DC power when 
round-robin scheduling is applied among multiple symmetric 



Eig. 5. The average output DC power when round-robin scheduling is 
performed among homogeneous (‘horn’) or heterogeneous (‘het’) UEs. 


or asymmetric UEs, as a function of the beacon transmit 
power. The transmit power (X-axis) and the average out¬ 
put DC received power (Y-axis) are shown in dBm, i.e. 
P = 101ogio(^o/10“^) where Pq is in the unit of watt. For 
illustration purposes, the distances of the symmetric UEs are 
set to 2m (i.e. at the inner boundary of the annular coverage 
as shown in Fig. [^, 10m and 30m (i.e. at the outer boundary 
of the annular coverage). It is observed from Fig. that the 
average output DC power decreases sharply with the increase 
in distance and that the simulation results are in perfect 
agreement with the numerical results computed by using ( [2Q| ). 

On the other hand, as shown in the dashed plot and with 
the plus/square marks, the heterogeneous scenario harvests 
the same average power as the symmetric UEs located at 
7.75m (computed by ( [3^ ) away from the power beacon, 
which coincides with the numerical results computed by ( [34| ). 
Consequently, for any particular coverage with inner and outer 
boundaries, the analytical expression shown in ( |^ (or ( [34| )) 
can be applied to efficiently predict the average output DC 
power in homogeneous (or heterogeneous) scenario, if round- 
robin scheduling is applied. Also, the effect of channel variety 
among asymmetric UEs can be exactly described by the 
distance scaling factor shown in ( [38] ). 


B. Opportunistic Scheduling 

With opportunistic scheduling performed among symmetric 
UEs, Fig. shows the average output DC power in linear 
scale, i.e. in the unit of mW (ImW = 10“^watt), where 
the beacon transmit power is set to 50dBm and the clustered 
UEs are located 10m away from the power beacon. It is 
clear that the simulation results of the average DC power 
agree exactly with the numerical results computed by ( [28] ). 
On the other hand, the scaling law given by ( [^ reflects 
the increasing trend of the DC power very well and, more 
importantly, all three curves demonstrate that the average DC 
power increases logarithmically with the number of UEs, as 
predicted by Theorem [^ Finally, we note that the small gap 
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Fig. 6. The average output DC power (in the unit of mW) when opportunistic 
scheduling is performed among homogeneous UEs (d = 10m), with the 
transmit power of the power beacon being 50dBm. 



Fig. 7. The average output DC power when opportunistic scheduling is 
performed among heterogeneous UEs, in the unit of mW. 


between the simulation and the numerical results of the scaling 
law, as shown in Fig.[^ is due to ignoring the Euler’s constant 
in ^2^ when ( [^ was derived. Also, for the special case of 
single UE, i.e. when N = 1 (corresponding to the most-left 
point of the lowest curve in Eig.[6|)^e numerical result of the 
scaling law was computed as perj^, since, in such a special 
case, the scaling law shown in ( [29] ) degenerates to zero. 

On the other hand, if opportunistic scheduling is performed 
among multiple asymmetric UEs, Eig. [^illustrates the average 
output DC power in the unit of mW, with respect to the number 
of UEs N. It is observed that, for either values of the beacon 
transmit power, 43dBm or 50dBm, the power scaling law 
computed by ( [47] ) perfectly reflects the increasing trend of the 
simulation results, i.e. linearly increasing with N. Moreover, 
when the transmit power is increased from 43dBm to 50dBm, 
if we show the output DC power in dBm, then it is easy to And 



Fig. 8. The power transfer efficiency. 


that the average output power corresponding to the case with 
the beacon transmit power set to P = 50dBm is exactly 7dB 
larger than that of the case with P = 43dBm, i.e. the same 
amount of increase as in the transmit power of the beacon. 
This is not strange since the output DC power at a particular 
receiver always increases linearly with the beacon transmit 
power, as shown in ( p^ . 

Compared with Eig. [^ it is evident that the curves in Eig. [7] 
have steeper slopes. This is because opportunistic schedul¬ 
ing yields larger diversity gain than round-robin scheduling, 
namely N versus In A^, as stated in Corollaries 1 and 2. 

C. Power Transfer Efficiency 

Based on the results shown in Eigs. and [7] the power 
transfer efficiency when opportunistic scheduling is applied 
among symmetric UEs is shown in the upper panel of Eig. 
and that among asymmetric UEs is displayed in the lower 
panel, in comparison with round-robin scheduling. The power 
transfer effici encies pertaining t o the round-robin mode are 
computed by Pjljg^/Perp, where is given by ( |20| ) in the 
symmetric case and by ( [34| ) in the asymmetric one. 

Erom the upper panel, it is observed that, the WPT effi¬ 
ciency with opportunistic scheduling increases logarithmically 
with the number of UEs while that of round-robin scheduling 
is constant at about 0.013%. On the other hand, the curves 
in the lower panel show that the WPT efficiency with oppor¬ 
tunistic scheduling increases linearly with the number of UEs 
while that of round-robin scheduling remains about 0.18%. 
In particular, when the number of UEs > 35, the WPT 
efficiency under opportunistic scheduling is higher than the 1% 
value which serves as a benchmark to the success of far-fleld 
WPT. Therefore, opportunistic scheduling is demonstrated to 
be a powerful technique to signiflcantly improve far-fleld WPT 
efficiency. Einally, it is worthwhile to point out that the ap¬ 
pearance of the WPT efficiency under round-robin scheduling 
looking higher than that under opportunistic scheduling when 
N is smaller than 5 (cf. lower panel of Pig. is due to the 
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fact that the average output DC power shown in ( [47] ) is an 
asymptotic result, which is not accurate when N is small. 

VII. Conclusion 

This paper provided a comprehensive study on the transmis¬ 
sion efficiency of far-field wireless power transfer (WPT), by 
taking into account both large-scale path loss and small-scale 
multi-path fading over bounded coverage area. The analyti¬ 
cally computed time-average output DC powers with respect to 
round-robin scheduling and opportunistic scheduling demon¬ 
strate that the latter scheduling strategy can improve far-field 
WPT efficiency linearly with the number of user equipments. 
This result provides a new approach to significantly improve 
the end-to-end power transfer efficiency of far-field WPT 
systems, apart from the traditional approach of increasing the 
transmit and receive antenna gains. 


Appendix A 
Proof of Lemma[T] 


we 


Before detailing the proof of Lemma in Section IV-B 
recall one of the seminal von Mises criteria | [2Q| , which will 
be used later in the form of Lemma lAll 


Lemma A1 (von Mises sufficient conditions for the domain of 
attraction of Gumbel distribution). Let F{x) be a distribution 
function, if there is a real number xi such that, for all xi < 
X < uj{F) where uj{F) = sup{x : F{x) < 1}, f{x) = 
F'{x) > 0 , and 


fix) 


, 

x^uj(F) 1 - F{x) 


= Co, 


(49) 


where cq G (0, oc) is a constant, then, F{x) is in the domain 
of attraction of the Gumbel distribution function Hs,o{x) = 
exp (—e“^). 


Applying Lemma 


A1 


to the CDF Fp 

X O 


t-i{t)ix) ([T^, it is 


clear that the parameters used in Lemma |A1| can be given 
by = 0 and uj{F) = oc. Then, when 0 < x < oc, the 
first-order derivative of ^^ists, i.e. /Pout-i(t)(^) 

in and is greater than zero. Hence, the first condition of 
the von Mises criterion is satisfied. Further, by virtue of the 
PDF in @ and the CDF in it is easy to show that 

/Pout-i(t) (x) 


lim 

x^oo 1 — Fid 

X o 


t-.(t)(x) 


cpa'h 


2 ’ 


(50) 


which implies that the second condition of the von 
Mises criterion is also satisfied. Consequently, by recall¬ 
ing Lemma A1 given by is in the 

domain of attraction of the Gumbel distribution function 
H^^q{x). That is, the limiting distribution function of the 
maximum among N output DC powers, i.e. Pout-?(^) = 
max{Pout-i(^),^out- 2 (t), • • • ,Pout-Ar(t)}, is given by the 
Gumbel distribution as shown in ( [25] ). Furthermore, the scaling 
factor ai is determined by Theorem 2.1.3] 


ai 


= inf 


1 — Fp 


,_,(t)(a^) < 2|. (51) 


Since the CDF given by is monotonically increasing, 
substituting ([T9J into •HD and performing some algebraic 


manipulations yields ( |^ . The positioning parameter hi is 
determined by pp. 104-105] 


/Pout-i(t)(®) 


(52) 


Substituting (p^-(p^ into ( [5^ results in ( [27] ), thus complet¬ 
ing the proof. 


Appendix B 
Derivation of Eq. 

We first derive the CDF and the PDF of df^{t). By 
definition, the CDF of df^{t) can be given by 

U-“(t)(^) = < a;} = 1 - , (53) 

where rex < x~i < Vnet, i-e. r~^ < x < r~^. Substituting 
into ( |5^ and performing some algebraic manipulations, 
we obtain 




^2 _ 2 
' net ' ex 


if J-ne? <X< 


(54) 


0, otherwise. 

Then, differentiating ( [54] ) with respect to x yields the PDF: 


r f \ ) a(r^ ) X - \ if “ < a; < 


(55) 

Now, we derive the CDF of the output DC power 
^out-i(^) = cl3df^{t)\hi{t)\‘^. By definition, we have 

= Pr {c/?d““(t)|/ii(t)|^ < x} 

OO 

= J Pr {cj3d~°‘{t)\hi{t)f < x\\hi{t)f = y} 

0 

OO 




i iy) dy 


OO 


c/3 'ex 


c/3 net 


Jf\hM\^iy)dy+ f\hi(t)\^iy)dy(56) 




r 5 
J 


net 'p^ 

^2 °!„2 f\h,(t)p{y)dy 

' net 'ex 




-L 


c/3 ’"net 1 t X \ 

-2 -^ ^ f\h,(t)\2iy)dy, (57) 

rLt - \cPy/ ' 


c/3 ex 


l2 


where ( [54] ) was inserted into ( [56] ) to reach ( [57] ). Subsequently, 
in light of ( [T^ , the first integral term in ( [57] ) is computed as 


h = 


' net ' ex 


/ ' c/3 ^net 

<Xh 


exp 


dy 


exp 


cjSal 
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- exp 


cj3al 


(58) 


Also, the second integral term in can be computed as 


h = 




- r 1 


L 


cl3 ^net / X \ “ 1 

\cl3y) ' 


:/3 ' ex 

X 

cpal 


exp dy 


a cfia'i 


r 1 + -, 


Q(’ c/3cr^ 


(59) 


Finally, inserting ([58|)-([5^ into •HD and performing some 
algebraic manipulations yields the desired 


Appendix C 
Proof of Lemma[2] 


Before detailing the proof of Lemma we reproduce 
another von Mises criterion different from Lemma A1 in Ap¬ 
pendix]^ as the following Lemma [CT] pT) Theorem 10.5.2]. 


Lemma Cl (von Mises sufficient conditions for the domain of 
attraction of Frechet distribution). Let F{x) be a distribution 
function with uj{F) = oo, if f{x) = F'{x) > 0 for all large 
X and for some ( G (0, (X)), 


lim 

x^oo 


xf{x) 

1 - F{x) 


= c, 


(60) 


then, F{x) is in the domain of attraction of the Frechet 
function Ffi^^(x) = exp (—, x > 0 . 


Recalling the fact that the first-order derivative of the 
complementary incomplete Gamma function is given by 
^r(a, 6 x) = —he~^^{hx)~^^^ and taking the derivative 
of ( [3^ with respect to x, yields the PDF of the output DC 
power. Then, substituting ( [^ and its PDF into ( [60| ), and 
performing some algebraic manipulations, it is easy to show 
that C = Therefore, the limiting distribution function of 
the maximum among N output DC powers, i.e. Pout- 2 (^) = 
max{Pout-i(^),^out- 2 (t), • • • ,^out-Ar(t)}, is given by the 
Frechet distribution as shown in ( [39] ). Furthermore, the scaling 
factor 62 in ( [^ is determined by Eq. (10.5.6)] 

62 = inf ja; ; 1 - . (61) 

Due to the high complexity of Ppout-i(t) (^) shown in ( [3^ , the 
exact solution to 62 is mathematically intractable. However, an 
explicit approximation of 62 can be derived as follows. 

In view of ( [3^ , we have at the top of the next page. 


1 “ PPout-iit)ix) 

= [’’net (“^2 x) - exp (-ri x)] 

' net ' ex 



X 


r 


2 

+ -,riX 
a 


-r 



exp{-rix) 
- »’lx) ri x 


2ri 


exp(-r 2 x) 


- xly,) r2 X 


(62) 


where the asymptotic series of the incomplete Gamma func¬ 
tion, i.e. r(a,x) = e~^ x^~^ [l + + 0 (^)], was em¬ 

ployed to get ( |62| ). Moreover, since r 2 ^ ri as defined after 
( [3^ , the second term on the right-hand side of ( |62| ) can be 
further ignored (this will be compensated later), yielding 

exp(—ri x) 


1 — Fp 

X o 


2 r 


V^nety riX 

Then, by recalling the fact that the CDF Pp^^^_.(^)(x) is 
monotonically increasing with respect to x, substituting ( [63] ) 
into and performing some algebraic manipulations yield 

V^net, 


rib 2 exp (r 162 ) 


(64) 


Since ^ f with e = 2.7183 • • • being the Euler’s 

yF net / ^ 

number, holds in practice, then applying the series expansion 
of the Lambert function Eq. (4.13.5)] to ( [^ results in 
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(65) 


Notice that ignoring the second term in ( [62] ) results in 
the value of 62 being underestimated, since the function 
exp(—x)/x involved in the first term of ( [ 6 ^ is monotonically 
decreasing with x > 0. To compensate this approximation er¬ 
ror, 62 in ( [^ is empirically enlarged (based on the observation 
from ( [ 35 ] )) so as to obtain the desired ( [40] ). The effectiveness 
of 62 in ( [40] ) is corroborated by ensuing simulation results in 
Section VLB Finally, it is noted that the value of 62 is not 
unique and that different choices lead to different convergence 
speeds of the limiting distribution function Section 2.10]. 
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